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Universal Amplitude Ratios in the 3D Ising Model 
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We present a high precision Monte Carlo study of various universal amplitude ratios of the three dimensional 
Ising spin model. Using state of the art simulation techniques we studied the model close to criticality in both 
phases. Great care was taken to control systematic errors due to finite size effects and correction to scaling terms. 
We obtain C+/C- = 4.75(3), f +, 2nd/ f -,2nd ~ 1.95(2) and u* — 14.3(1). Our results are compatible with those 
obtained by field theoretic methods applied to the 0* theory and high and low temperature series expansions of 
the Ising model. 



1. Introduction 

In the neighbourhood of a second order phase 
transition various quantities display a singular be- 
haviour. In this limit microscopic features be- 
come irrelevant and models which differ at the 
microscopic level may share the same singular be- 
haviour. This is the basis of universality. Dif- 
ferent models belonging to the same universality 
class share the same critical indices. However uni- 
versality has much stronger implications and it is 
possible to show that for a given universality class 
the values of particular critical-point amplitude 
combinations are unique |§. In the following we 
shall be interested in the universality class of the 
3 dimensional Ising model which has several in- 
teresting experimental realizations, ranging from 
the binary mixtures to the liquid vapour transi- 
tions. The action of the Ising model is given by 



(1) 



<n,m> 



where s„ G { — 1, 1} is the field variable and 
< n,m > is a pair of nearest neighbour sites on 
the lattice. f3 = is the coupling. The reduced 



temperature t can be written as t 



where 



/3c is the critical coupling. In the following we con- 
sider lattices of volume and periodic boundary 
conditions. 



We simulated the model in the low as well as in 
the high temperature phase. The simulations in 
the high temperature phase have been performed 
using Wolff's single cluster algorithm. Here the 
improved estimators give a great boost to the ac- 
curacy of the results. However in the low tem- 
perature phase, due to the finite magnetisation, 
the improved estimators of the cluster-algorithm 
are of little help. Hence wc simulated here with 
a multi-spin coded demon algorithm. 

We measured the magnetisation, the magnetic 
susceptibility x s-nd the exponential ^ as well as 
the second moment correlation length ^2nd- In or- 
der to obtain a finite result for the magnetisation 
in the low temperature phase on finite lattices at 
vanishing external field we measured 
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The magnetisation defined in this way was also 
used for the subtraction in connected correlation 
functions in the low temperature phase. 

In the neighbourhood of the critical point one 
expects a singular behaviour of the form m ^ 
B tf, X - C±\t\-\ C - /il^r" and ~ 
f±,2nd\t\~'^ , where -I- indicates the high temper- 
ature phase t > and — the low temperature 
phase t < 0. 

Ratios and more general combinations of sin- 
gular quantities where the singularities cancel are 
expected to be universal. In our study we com- 
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puted f+^2nd/f^.2nd, 3 2„,B2 

and C+/ f^l 2nd^^ ■ The finiteness of the last two 
quantities rehes on the scahng (hyper-scahng) re- 
lations a + 2/3 + 7 = 2 and dv — 2 — a. 



2. numerical results 

First we simulated the model in the low tem- 
perature phase at /? = 0.2391, 0.23142, 0.2275, 
0.2260,0.2240 and 0.22311 on lattices of size L = 
30 up to L = 120. From finite size tests at the 
largest /3 we concluded that L > 20^ is needed to 
be save of finite size effects. The total number of 
measurements was about 3 x 10® for all /3's. We 
carefully compared our results with IDA's of low 
temperature series expansions and Monte Carlo 
results given in the literature. In general we found 
good agreement. For details see ref. po[| . 

Next we simulated the model in the high tem- 
perature phase at /3 = 0.20421, 0.21189, 0.21581, 
0.21731, 0.21931 and 0.22020. In the high tem- 
perature phase only L > 6^ is required to be suf- 
ficiently close to the thermodynamic limit. The 
/3- values were chosen such that they pair up with 
a /3 in the low temperature phase such that 

Plow + Phigh — 2/3c- 

This choice of /3-values allows to compute for 
example approximations of C+ / C_ by Ty^ (t) — 
-^^[j with t > 0. In this way we need not to ex- 
tract the C± themselves and the critical exponent 
7 does not enter into the calculation. The critical 
limit is then obtained from a fit to the ansatz 



T^{t) = C+/C- + ao r 



fli t 



(3) 



where we have included the leading corrections 
to scaling. The best know numerical estimate of 
the correction exponent is 6 = 0.51(3). The re- 
sults for the other amplitude ratios and combi- 
nations are extracted analogously. In table 1 we 
have summarised the approximations of C+/C_, 
f+,2nd/f^,2nd and u* = 3C_//i 2nd-S^ ^nite 
t. The last line gives our result for the extrapo- 
lation to the critical limit. In addition we have 
C+//+,2nd^^ = 3.05(5) The largest reduced tem- 
perature t is excluded from the fit in all cases. 
The errors quoted for the critical limit include 
errors induced by the uncertainty of Pc and 9. 
The x^/d-o.f. was of order 1 for all quantities. 



Table 1 

The various ratios as functions of the reduced 
temperature t The last hne gives the extrapola- 
tion to the critical limit 

t Fv Ft u 



0.0787 
0.0441 
0.0264 
0.0196 
0.0106 
0.0066 



6.044(5) 
5.546(5) 
5.283(3) 
5.182(5) 
5.027(6) 
4.947(6) 



1.902(2) 
1.920(3) 
1.932(3) 
1.939(3) 
1.942(3) 
1.948(3) 



15.00(6) 
14.75(5) 
14.66(6) 
14.64(5) 
14.47(6) 
14.51(7) 







4.75(3) 1.95(2) 14.3(1) 



Table 2 

Results for the amplitude ratios reported in the 
literature. They are obtained by e-expansion [2], 
perturbation theory in 3 dimensions [2,3], high 
and low temperature series expansions of the 
Ising model [4] and Monte Carlo simulations [5]. 



ref. 



f- 



[2\ 4.70(10) 

[21 4.82(10) 

[31 4.72(17) 2.013(28) 14.4(2) 

[i| 4.95(15) 1.96(1) 14.8(1.0) 

[5| 5.18(33) 2.06(1) 17.1(1.9) 



3. Comparison with theoretical results and 
experimental data 

In table 2 we have summarised the most re- 
cent results for amplitude ratios and combina- 
tions for the 3D Ising universality class obtained 
from e-expansion |^] , perturbation theory in 3 di- 
mensions [HJsl, high and low temperature series 
expansions of the Ising model and a Monte 
Carlo study of the Ising model . The results 
for C+/C- of refs. [g-~|§[ are in good agreement 
with our estimate. Our result is however more 
accurate. For f2nd.+/ f2nd,- our result is compat- 
ible with that of Q[ , the result from perturbation 
theory in 3 dimensions [^ is larger then our esti- 
mate, however still within a 2a deviation. There 
is however a clear discrapency with the Monte 
Carlo result of ref. a . Their value is much larger 
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then our one. For u* and -73 — ^-^2 our result is 

consistent with the other methods. In both cases 
we could reduce the errors considerably. 

The experimental data reported in tab. 3 refer 
to the three most important experimental realiza- 
tions of the Ising universality class, namely binary 
mixtures (bm), liquid- vapour transitions (Ivt) 
and uniaxial anti- ferromagnetic systems (af). It 
is important to notice that these realizations are 
not on the same ground. Antiferromagnetic sys- 
tems are particularly apt to measure the C+/C_ 
and f +. 2nd / f -,2nd ratios, while for the liquid 
vapour transitions the Fc = Rc/ combination 
is more easily accessible. Finally, in the case of 
binary mixture all the three ratios can be rather 
easily evaluated. The common attitude is to as- 
sume that the above systematic errors are ran- 
domly distributed and to take the weighted mean 
of the various experimental results. The numbers 
are taken from ref . ^ , where more details on the 
experiments and the averaging procedure can be 
found. 



Table 3 

Experimental estimates for some amplitude ra- 
tios, (bm) refers to binary mixtures, (Ivt) to 
liquid-vapour transitions and (af) to uniaxial 
anti- ferromagnetic systems, (all of them) gives 
an weighted average. 



exp. setup 


C- 


f+,2nd 


0+ 

n.2^dB'' 


(bm) 


4.4(4) 


1.93(7) 


3.01(50) 


(Ivt) 


4.9(2) 




2.83(31) 


(af) 


5.1(6) 


1.92(15) 




(all of them) 


4.86(46) 


1.93(12) 


2.93(41) 



In the case of the f +, 2nd / f -,2nd ratio (for 
which, as we have seen, some of the present the- 
oretical or Monte Carlo estimates disagree) we 
have listed, for a more detailed comparison, all 
the available experimental data in tab. 17. In this 
table we denote with "N-H" the nitrobenzene - n- 
hexane binary mixture, and with "I-W" the one 
obtained by mixing isobutyric acid and water. A 
much more detailed account of the various exper- 



imental estimates can be found in |l 



Table 4 

Experimental estimates for the {^'^"'^ ratio 



Ref. 


year 


exp. setup 


/+,2Tld 

f-.2nd 


1 
1 
1 
1 


6 
7 

S 
9 


1 
1 
1 
1 


1972 
1980 
1983 
1986 


(af),FeF2 
(af),CoF2 
(bm),N-H 
(bm),FW 


2.06(20) 
1.93(10) 
1.9(2) 
2.0(4) 



A detailed account of our work can be found in 
ref 0. 
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